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■ To describe high energy collisions one widely uses thermodynamical methods and con- 

■ cepts which follow the classical Boltzmann-Gibbs (BG) approach. In many cases, however, 
^p H ' either some deviations from the expected behaviour are observed experimentally or it is 

known that the conditions necessary for BG to apply are satisfied only approximately. 

i-C ■ In other branches of physics where such situations are ubiquitous, the popular remedy 

is to resort, instead, to the so called nonextensive statistics, the most popular example 

^ ■ of which is Tsallis statistics. We shall provide here an overview of possible imprints of 

non-extensitivity existing both in high energy cosmic ray physics and in multiparticle 
production processes in hadronic collisions, in particular in heavy ion collisions. Some 
novel proposition for the interpretation of the nonextensitivity parameter q present in 
such circumstances will be discussed in more detail. 
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1 Introduction 



It has been realized for some time already that in many places of science there are phenom- 
ena which clearly indicate the existence of some degree of nonextensitivity (understood in the 
thermodynamical sense). They include all situations characterized by long-range interactions, 
long-range microscopic memory and space-time (and phase space as well) (mult i) fractal struc- 
ture of the process. Such anomalous (from the point of view of the standard thermodynamics 
and Boltzmann-Gibbs statistics) systems are found to be best described in terms of generalised, 
nonextensive thermostatistics, the most popular and explored example of which is the so called 
Tsallis statistics which is characterised by the nonextensitivity parameter q [|1|] . 

To make our presentation self-contained we shall first provide the basic formulas (refering to 
cf. [|IJ and references therein for a thorough discussion of all possible aspects of nonextensivity) . 
Everything is based on a generalized entropic form depending on a single parameter (entropic 
index) q in such a way that for q — > 1 it gives the normal BG entropy: 

s q = ~yzt f 1 _ 2?? ) ^ s BG = - Y^Pi^Pi- (!) 

The S q is nonextensive in the sense that 

S q (A + B) = S q {A) + S q (B) +{l-q) S q {A)S q (B), (2) 

where A and B are two independent systems in the usual sense, i.e., Pij{A + B) = pi(A)pj(B). 
In this sense the entropic index q is also a measure of the nonextensivity in the system. Using 
the usual procedure of information theory when looking for the most probable and least biased 
(normalized) probability distribution of some events x subjected to a single constraint in the 
form of unnormalized g-expectation value < A > q = J dxA(x)p(x) q [fj], one gets immediately, 
from the maximization of the entropy S q the known expression 

p q {x) = i-[l - (l-q)aA(x)}^ ^ i exp (-aA(x)) . (3) 

The Lagrange multiplier a can be expressed in terms of the < A > q from the imposed constraint 
and Z q from the normalization condition. Notice that whereas in the extensive, i.e., q = 1 case, 
all values of x G (0, oo) are admissible, for nonextensive case of q ^ 1 we have restrictions so 
that [1 — (1 — q)aA(x)] is positive. 

Our presentation will be devoted to a very limited subject of high energy collisions. We 
shall attempt to overview the probable imprints of nonextensivity in high energy collisions 
where under this term we understand both multiparticle production processes taking place in 
cosmic ray experiments and those observed in accelerator experiments. 

In cosmic ray experiments one encounters routinely a cascade processes, both in the atmo- 
sphere and in emulsion chambers serving as detectors. The former originate during the passage 
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of the primordial cosmic rays through the atmosphere with subsequent collisions and dissipation 
of energy (they are known as the so called Extensive Air Showers) - this is typical stochastic 
process not very much sensitive to details of the elementary interactions ||. The latter are 
connected with the actual detection process taking place in special emulsion chambers exposed 
in many places of the Earth, usually at high altitudes. Some of their characteristics do depend 
on details of the interaction process. Because in both cases one encounters formulae of the type 
(|3]), they are a priori sensitive to the possible nonextensivity of such processes. 

The accelerator high energy collisions are usually connected with production of large num- 
ber of secondaries (mostly ir and K mesons) . The strong interactions involved here make their 
detail descriprion from first principles impossible and one is forced to turn to phenomenological 
models of various kinds. The most economical (as far as the number of parameters is concerned) 
are thermodynamical and statistical models which have been, in fact, in use since the beginning 
of this field of research. It should be stressed that the very first successful phenomenological 
model of multiparticle production, the so called Landau Hydrodynamical Model, was proposed 
already in the pre-accelerator era of 1953 and used in the analysis of specific multiparticle data 
taken from cosmic ray interactions . This model paved a way to more sophisticated statistical 
models now in use, mainly in high energy heavy ion collisions which are believed to lead to the 
production of new state of matter, the Quark Gluon Plasma (QGP) Q. 

To recapitulate: in both cosmic ray and accelerator experiments one observes multiparticle 
production. However, in the cosmic ray case one also deals with propagation of the secondaries 
originating in multiple production processes at a given point and with their subsequent sec- 
ondary interactions (also of multiparticle type), i.e., with a full fledged cascade process ||. 

In Section 2 the possible traces of nonextensivity apparently seen in some cosmic ray experi- 
ments with emulsion chambers will be presented along with similar effects seen in elementary 
e + e~ H and hadronic collisions and in the collisions of nuclei in accelerators || [10| . In all of 
them distributions given by Eq. (||) are observed. Section 3 will be devoted to disscusion on 
possible hints emerging from the occupation number distributions < n > q ||, [10| . In Section 
4 we shall propose a novel interpretation of the entropic index q Jll| encountered in Sec. 2. 
The last Section contains final remarks together with a list of other possible hints of nonexten- 
sivity not discussed in detail here (including topics from QGP physics [[L2[] and nonextensivity 
manifested in the statistics of the quantum states produced in the scattering process such as 
7t— nucleon or n— nucleus scatterings fT3|| ). 
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2 Traces of nonextensivity in nonexponential distribu- 
tions 



Our encounter with the notion of nonextensivity started when we realised that one of our previ- 
ous results concerning the occurence of the so called long flying component in the propagation of 
the initial flux of incoming cosmic ray particles (mostly nucleons) |14| can be interpreted as yet 



another manifestation of the Levy distribution (pi) M . To be more specific let us briefly summa- 



rize the result of ||14|| . We have analysed, there, the distribution of cascade starting points in the 
extrathick lead chamber of the Pamir experiment. The corresponding data points for the num- 
ber of cascades originating at depth T (measured in cascade units, 1 c.u. = 6.4g/cm 2 = 0.56cm) 
are shown in Fig. 1. Whereas at small depths (up to ~ 60 cm of lead) we observe the usual 
absorption of hadrons as given by the simple exponential formula 

—— = const • exp I — — ) , (4) 



dT V A. 

at biggest thickness there is, against all expectations, noticeably excess of experimental points 
above the simple extrapolation of small-depth data. The observed discrepancy means that orig- 
inal hadrons tend to fly longer without interaction (that is why a term long flying component 
is coined for this type of phenomenon). 



In |14| we have argued that the observed effect can be just another manifestation of the 
fluctuation of the corresponding hadronic cross section a = Am^ j (where A denotes the mass 
number of the target and is the mass of the nucleon, such a possibility is widely discussed 



in the literature and observed in diffraction dissociation experiments on accelerators, cf. [14 
for details and references). It turned out that fluctuations of this cross section (i.e., in effect, 
fluctuations of the quantity 1/A) with relative variance 

u =U-LL>o.2 (5) 
allow to describe the observed effect. 

It turns out |7[ that the same data can be fitted by the nonextensive formula 

(6) 



dN 

—— = const 
dT 



1 -(!-«>! 



i 
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with parameter q = 1.3 (in both cases A = 18.85 ± 0.66 in c.u. defined above), cf. Fig. 1 [15 



Similar example is also known in heavy ion collisions ]T0[ (see also 0). It turns out that dis 



tributions of transverse momenta of secondaries produced in nuclear collisions at high energies 
(transverse with respect to the collision axis given by the direction of the colliding objects in 
the center of mass frame) dN(pT)/dpT are best described by a slightly nonexponential function 
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of the type 



dN(p 7 
dp T 



const 



(1-9) 



m 2 + p\ 
~~kT 



i 
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const • exp 



m? + p\ 



kT 



• (7) 



Here m is the mass of produced particle, k is the Boltzmann constant (which we shall, in what 
follows, put equal unity) and T is, for the q = 1 case, the temperature of the reaction consid- 
ered (or, rather, the temperature of the hadronic system produced). In fact, precisely from 
such exponential fits to transverse masses my = 



'm 2 + pj, one infers information about the 
temperature T. Therefore, any deviation from the exponential behaviour of such distributions 
are always under detailed scrutiny in which one is searching for the possible causes. In || it 
was suggested that the extreme conditions of high density and temperature occuring in ultra- 
relativistic heavy ion collisions can lead to memory effects and long-range colour interactions 



and to the presence of non-Markovian processes in the corresponding kinetic equations [16 



As it is seen in Fig. 2 one indeed finds [KJ a small deviation from the exponential 



behaviour (on the level q = 1.015). As we shall demonstrate in Section 4 it can, however, lead 
to quite dramatic effects. It was also shown in |§ that to first order in \q — 1| the generalized 
slope becomes the quantity 

T q = T + {q- l)m T . (8) 

with T being temperature of a purely thermal source. This should be contrasted with the 
empirical relation for the slope parameter T, from which the freeze-out temperature (at which 
hadrons are created from the QGP) Tf is then deduced, 

T = Tf + m(v ± } 2 . (9) 

The (v±) is a fit parameter usually identified with the average collective (transverse) flow veloc- 
ity of the hadrons being produced. In @ one has, instead, a purely thermal source experiencing 
a kind of blue shift at high itlt (actually increasing with The nonextensivity parameter q 
accounts here for all possibilities one can find in [|16| and could, therefore, be regarded as a new 
way of presenting experimental results with q ^ 1 signaling that there is something going on in 
the collision that prevents it to from being exactly thermal-like in the ordinary sense mentioned 
above. 



We shall proceed, now, to two other examples of possible nonexponential distributions. First 
is the attempt || to fit the energy spectra in both the longitudinal and transverse momenta 
of particles produced in the e + e~ annihilation processes at high energies. Those are, contrary 
to the previous example, the most elementary high energy multiparticle production processes. 
The initial e + e~ pair annihilates to a virtual photon which subsequently gives rise to a (highly 
excited) quark-antiquark pair. They in turn develop a complex hadronization process related 
to the long-distance (strong coupling) regime of Quantum Chromodynamic (QCD). Usually 
being described in terms of the so called string model ]17| it admits also, for low energies, a 
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kind of thermodynamical equilibrium approach ||. However, it turns out that when going to 
higher energies one cannot keep the temperature T inferred from the px distributions constant, 
invalidating therefore the whole concept. On the other hand, using instead the nonextensive 
power-like (g-dependent) distribution, one can write the following transverse momentum dis- 
tribution 101: 



1 da 
a dpx 



const ■ p T / dp i 
J o 



pj + m z + 

1 - (l-g)^^— f 

-in 
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(where p^ is longitudinal momentum of secondary particle of mass m and q is the entropic 
index). Keeping the temperature To essentially constant and changing only q one can now fit 
data extremely well, cf. Fig. 3 [BJ. In this example the q ^ 1 is then regarded as a manifesta- 
tion of nonextensivity arising in the hadronization process in which quarks and gluons combine 
together forming hadrons, a process which involves long range correlations in the phase space. 
Actually, this observation has general validity and applies to all production processes discussed 
here as well. It applies, most probably, also to a pure elastic quantum scattering processes 
discussed in fll3||. 



The final example from this category deals with the most probable rapidity distributions for 
produced secondaries. Suppose that we have an excited object of mass M which hadronizes into 
N secondaries of transverse mass mx each. For simplicity we shall consider only one-dimensional 
hadronization in which transverse momenta are hidden in the mean (px) parameter, which is 
kept constant and enters In this case the observable distribution of interest is f(y) = jf^r~> 
i.e., a distribution in rapidity defined as y = \ ha[(E + pl)/{E — Pl)}- Some time ago it was 
shown, using the maximization of the information entropy in the Shannon form (i.e., in fact 
the BG one), that knowing only the energy M and multiplicity N one has to expect that |L8 



f(y) = z(Af N) exp [—/3(M, N) ■ m T coshy] , (11) 

where Z(M,N) comes from the normalization of f(y) and the Lagrange multiplier (3(M,N) 
from the energy conservation constraint. Of special interest to us is the fact that for some 
values of the mean energy M/N (i.e., for some values of N for a given mass M) j3 can be zero 
and even negative. Actually it can be shown that (3 > only if 

M 

N > N ~ 2 In . 12 

m T 

This statement invalidates the widely assumed, on different occassions, the so-called Feynman 
scaling hypothesis, namely, that in high energy collisions one should expect f(y) ~ const (cf. 
[]l8l for more information). Such hypothesis is clearly incompatible with the experimentally 
observed fact that the mean multiplicity < N > grows with energy faster than In M, more like 



< iV >~ M 0A °' 5 ). However, using Tsallis g-entropy, instead of BG one obtains |L9 



fq{y) = Z (M, N) [1 " (i-^/W^-^TCOsh^ (13) 



6 



with P q (M, N)>0 for 

/ M \ q 

N >N q0 ~2 In . (14) 

V m T J 

It is, then, obvious that for q > 1 one can indeed accommodate, at a given rapidity interval, 
more particles with f(y) = f q (y) being constant. In a sense one can think of a kind of Feynman 
q-scaling here (which would generalize the usual one and describe, in terms of the parameter q, 
what is usually called a violation of the Feynman scaling hypothesis |2(|). In Fig. 4 this effect 
is clearly demonstrated for M = 100 GeV (and = 0.4 GeV). While in Fig. 4a N = 10 is 
chosen in such a way as to have /3(100, 10) ~ (notice that in this case /3 g= o.7 > whereas 
/3 g= i.3 < 0), in Fig. 4b one has, instead, /3 g= i. 3 (100, 20) = (and remaining [3 q= i > and 
P q =o.7 > 0). Fig 4c shows, for comparison, the case with N > N Q for which all (3 q > 1. 



3 Traces of nonextensivity in the mean occupation num- 
bers n q 

Another place where nonextensivity enters in a natural way is the mean occupation numbers 
generalizing the Bose-Einstein or Fermi-Dirac ones to a q ^ 1 case. Whereas the single particle 
distribution function is obtained in the usual procedure of maximizing the Tsallis entropy under 
the constraints of given average internal energy and number of particles, the mean occupation 
numbers (n) q are not available in analytical formula for any q. Only in the dilute gas approx- 
imation and for small deviations of q from unity can one express them in a simple analytical 



form 21 



(n) q = {[1 + (g - l)fi{E - /i)] 1 ^ ± I}"' , (15) 

where (3 = 1/kT, [i is the chemical potential and the +/— sign applies to fermions/bosons. 
Notice that in the limit q — > 1 (extensive statistics) one recovers the conventional Fermi- 
Dirac and Bose-Einstein distributions. What will interest us here are the generalized particle 
fluctuations, 

(^^^-Tu^m^7) ilTin) ' h (16) 



where E = y/m 2 + p 2 . That is because, so far, these formulas have been applied to study the 



fluctuation pattern expected in heavy ion collisions |9|, |10| in measurements performed on an 
event-by-event basis. Notice that the denominator occuring in flT6|) modulates in a novel and 
specific way the usual pattern of fluctuations for the q = 1 case 0. 

Event-by-event fluctuations can be used as a valuable source of information on the dynamics 
of heavy-ion collisions. However, there is the problem of how to disentangle the dynamical 
fluctuations of interest from the trivial geometrical ones due to the impact parameter variation 
(resulting in the different number of nucleons participating in a given event). To solve this 
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problem it was proposed in 
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to use the following measure of fluctuations or correlations: 



N 



where Z = 2J ^i- (17) 



i=l 



Here where x denotes the mean value of the observable x calculated for all particles 

from all events (the so called inclusive mean) and N is the number of particles analysed in the 
event. In (|17|) (N) and (Z 2 ) are averages of event-by-event observables over all events whereas 
the last term is the square root of the second moment of the inclusive z distribution. $ equals 
zero when the correlations are entirely absent. On the other hand, it is constructed in such a 
way as to be exactly the same for nucleon-nucleon and nucleus-nucleus collisions if the latter is 
a simple superposition of the former. The ^-measure has been successfully applied to the ex- 
perimental data (cf. |p3[| ) and the fluctuations of transverse momentum px, which are observed 
in nucleon-nucleon collisions, have been found to be significantly reduced in the central Pb-Pb 
collisions at 158 GeV per nucleon. Its nonextensive extension has been analysed in || where 
it was found that the usual correlations are increased for q < 1 and decreased for q > 1. This 
time the measure $ can became positive or negative for both fermions and bosons depending 
on the value of the entropic index q. Actually in || it was found that for T = 140 MeV and 
fi = the ^-measure vanishes for q = 1.015. 

The clear prediction of M, which can consist a subject for experimental verification, is that 
the pr-dependence of partial contributions to $ PT should become negative for px > 0.5 GeV. 
Notice that in this regime the ^-measure should already be free from contaminations from res- 
onance decays and, therefore, the experimental confirmation of this prediction would provide a 
strong signal for the nonextensivity present in relativistic heavy ion collisions. 

The ^-measure is applicable not only to fluctuations of kinematical quantities such as pt 
but also to the azimuthal p4| and chemical fluctuations as well. The latter were analysed in [25] 



for the normal statistics and in JT0[ for the nonextensive one. The representive sample of results 



is shown in Fig. 5 for q = 1.015 mentioned above. For simplicity we have restricted ourselves, 



here, only to comparison with results of |25[ without resonances ||26| . Actually, one expects 
that for a given q fluctuations should grow with the mass of detected particle - this observa- 
tion provides yet another possibility for experimental verification of the nonextensivity concept. 



4 Nonextensivity parameter q as a measure of fluctua- 
tions 

The general picture emerging from the previous discussion is that one can account very eco- 
nomically (by introducing only one new parameter q) and adequately (by using nonextensive 
formulas emerging from Tsallis statistics with entropic index q) for a number of observations 
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deviating from the normal BG approach. The question of the possible meaning of entropic in- 
dex in these cases is therefore very natural. For the cases discussed in Section 2 we would like to 
propose that q is connected with fluctuations present in the system under investigation. Notice 
that common feature of the first two examples in Section 2 is that they both are described by 
the powerlike distribution of the type 



L q (e) 



a, 



A 



1-9 



As mentioned before, the cosmic ray example was originally explained [[bf] by the apparent 
fluctuation of the mean free path parameter A in the corresponding exponential formula 
It is then natural to expect that these fluctuations (which were so far described in |14| only 
numerically by means of Monte Carlo simulations) should be formulated in such a way as to 
result in eq. ([18]) with a parameter q. The same should be also true for the heavy ion collision 
example. Actually, this example is even more important and interesting because of the long and 
still vivid discussion on the possible dynamics of temperature fluctuations and 
because of its connection with the problem of QGP production in heavy ion collisions f28|, |30 



We shall, therefore, treat both representing the same class of fluctuation phenomena 

and claim that the parameter q is a measure of fluctuations present in the Levy-type distribu- 



tions ( p~8|) describing the particular process under consideration. 



To demonstrate this conjecture let us analyse the influence of fluctuations of the parameter 
1/A present in the exponential formula L q= i(e) ~ exp(— e/X). Our aim will be to deduce the 
form of a function /(1/A) which transforms the exponential distribution to a power-like Levy 
distribution flTB| ) and which describes fluctuations about the mean value 1/Ao- Although in 
both examples considered above the data preferred q > 1, we shall discuss q < 1 case as well. 
In the q > 1 case, where e G (0, oo), one has, 



Lq>l{£] A ) 



a, 



e 1 > ~ Q 

1 + T~ 

Aq ci, 



a, 



(19) 



where a = Writing the following representation of the Euler gamma function |31] 



1 + 



e 1 
An a 



T(a) 



drj rj 



Q-l 







e 




exp 


-v (i - 




a). 







(20) 



and changing variables under the integral to rj 
by the following gamma distribution: 



a 



An 

A ' 



one obtains eq. (O) with /(1/A) given 



f q >i (I 



fa 



1 1 

A' Ao" 



r(a) VA 



a-l 



cxp 



with fi = a\o and with mean value and variation in the form: 



— and 

Ao 




(21) 



a An 



(22) 
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Notice that, with increasing a the variance (|22|) decreases and asymptotically (for a — > oo, i.e, 
for q — > 1) the gamma distribution ([2~l"l) becomes a delta function, / (?>1 (1/A) = 5(X — A ). The 
relative variance for this distribution is given by 



UJ 



ay 



- = Q 

a 



1. 



(23) 



For the q < 1 case e is limited to e G [0, A /(l — q)\- Proceeding in the same way as 
before, i.e., making use of the following representation of the Euler gamma function (where 



a 



-a 



1 N 

l-q' 



[l \ ' 


a' 


( a'\ \ 


= 1 


QJ'Aq. 




\a'X -e) 


r(a') 



c??7 vj 



a'-l 



exp 



-77 1 



a'X — e 



(24) 



and changing variables under the integral to r] = "' A ° £ , we obtain L g<1 (£; Aq) in the form of 
eq. ( |I9D but with a — > —a' and with the respective /(1/A) = f q< i(l/X) given now by the same 



gamma distribution as in but this time with a — > a' and /i = = a'A — e. Contrary 
to the q > 1 case, this time the fluctuations depend on the value of the variable in question, 
i.e., the mean value and variance are now both e-dependent: 



1 



An 



and 




1 



a' 



(Ac 



(25) 



However, the relative variance 



— = i 

a' 



(26) 



remains e-independent and depends only on the parameter q. As above the resulting gamma 
distribution becomes a delta function, / 9<1 (1/A) = S(X — A ), for a' — > oo, i.e., for q — > 1. 

This completes the proof of our conjecture. The nonextensivity parameter q in the L q (e) dis- 
tributions can, indeed, be expressed by the relative variance uj of fluctuations of the parameter 
1/A in the distribution L q=1 (e): 

q = 1 ± uj (27) 

for the q > 1 (+) and q < 1 (— ) cases. 



Concerning transverse momentum distributions in heavy ion collisions, cIN(pt)/cIpt, it is 
interesting to notice that the relatively small value q ~ 1.015 of the nonextensive parameter 
obtained there ||, |10j, if interpreted in the same spirit as above, indicates that rather large 
relative fluctuations of temperature, of the order of AT/T ~ 0.12, exist in nuclear collisions. 
It could mean therefore that we are dealing here with some fluctuations existing in small parts 



of the system in respect to the whole system (according to interpretation of [p7|) rather than 
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with fluctuations of the event-by-event type in which, for large multiplicity N, fluctuations 



AT/T = 0.06/viV should be negligibly small f30|. This controversy could be, in principle, 
settled by detailed analyses of the event-by-event type. Already at present energies and nuclear 
targets (and the more so at the new accelerators for heavy ions like RHIC at Brookhaven, now 
commisioned, and LHC at CERN scheduled to be operational in the year 2006) one should 
be able to check whether the power-like pr distribution dN(pT)/dpx occurs already at every 
event or only after averaging over all events. In the former case we would have a clear signal of 
thermal fluctuations of the type mentioned above. In the latter case one would have for each 
event a fixed T value which would fluctuate from one event to another (most probably because 
different initial conditions are encountered in a given event). 

The proposed interpretation of q leads immediately to the next question: why and under 
what circumstances is it the gamma distribution that describes fluctuations of the parameter 
A? To address it let us write the usual Langevin equation for the stochastic variable A |j32| : 

nr + \- + £(*)! A = = const > 0. (28) 
at [t J 

with damping constant r and source term 0. This term will be different for the two cases 
considered, namely: 

1 ( 8 \ i j. ± Xo 



V = Vq<i = - Xo 7 whereas = <j) q>1 = — . (29) 

TV Or/ T 

For stochastic processes defined by the white gaussian noise form of £(£) one obtains the 
following Fokker-Plank equation for the distribution function of the variable A [[34] 



where the intensity coefficients K\ 2 are defined by eq.fl2"g|) and are equal to (cf., for example, 



35 



K\{\) = (p - - + D X and K 2 (\) = 2D\ 2 . (31) 

T 

From it we get the following expression for the distribution function of the variable A: 

(32) 



2 l X M K M 



/(A) = xfw exp L - a 

which is, indeed, a gamma distribution in variable 1/A, 

/(A) ^"©°" exp (-x)' (33) 

with the constant c defined by the normalization condition, / °° d(l / X) f (1 / X) = 1 and depending 
on two parameters: 

ri £ ) = < ^y- and aq = Vd' ( 34 ) 

with (j) q = g >i, g <i and a q = (a, a') for, respectively, q > 1 and q < 1. This means that we 
have obtained eq. (|2?D with uj = and, therefore, the parameter of nonextensivity q is given 
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by the parameter D and by the damping constant r describing the white noise. 



The above discussion rests on the stochastic equation (|28|). To comment on its possible 
origin let us turn once more to fluctuations of temperature |27|, [3(| discussed before, 

i.e., to A = T. Suppose that we have a thermodynamic system, in a small (mentally separated) 
part of which the temperature fluctuates with AT ~ T. Let X(t) describe stochastic changes 
of the temperature in time. If the mean temperature of the system is (T) = To then, as result 
of fluctuations in some small selected region, the actual temperature equals T' = T — r£(t)T. 
The inevitable exchange of heat between this selected region and the rest of the system leads 
to the equilibration of the temperature and this process is described by the following equation 

n 

^ - ~ (r - t) + n q = o (35) 

which is, indeed, of the type of eq. ( P8| ) (here Q q< i = -^r and fl q> \ = 0). 

In this way we have recovered eq. (^) and clearly demonstrated the plausibility of our 
proposition. Notice the presence of the internal heat source in the above equation in the q < 1 
case. It has a sense of dissipative transfer of energy from the region where (due to fluctuations) 
the temperature T is higher. It could be any kind of convection-type flow of energy; for exam- 
ple, it could be connected with emission of particles from that region. The heat release given 
by ej (to;') depends on e (but it is only a part of e that is released). In the case of such energy 
release (connected with emission of particles) there is additional cooling of the whole system. 
If this process is sufficiently fast, it could happen that there is no way to reach a stationary 
distribution of temperature (because the transfer of heat from the outside can be not sufficient 
for the development of the state of equilibrium). On the other hand (albeit this is not our case 
here) for the reverse process we could face the "heat explosion" situation (which could happen 
if the velocity of the exothermic burning reaction grows sufficiently fast; in this case because 
of nonexistence of stationary distribution we have fast nonstationary heating of the substance 
and acceleration of the respective reaction). 

It should be noticed that in the case of q < 1 the temperature does not reach stationary 
state because, cf. Eq. (|2lS|), (1/T) = 1/(T — e/a'), whereas for q > 1 we had < 1/T >= 1/T . 
As a consequence the corresponding Levy distributions are defined only for e G (0,T o a') be- 
cause for e —> T a', < T >— > 0. Such asymptotic (i.e., for t/r — > oo) cooling of the system 
(T — > 0) can be also deduced form Eq. ( |55| ) for e — > T a'. 



Our explanation, being tied to specific examples (especially to the example of the temper- 
ature fluctuations) differs from other works in which L q ^i(e) is shown to be connected with 
L g= i(e) by the so called Hilhorst integral formula (the trace of which is our eq. fl2~0D) |3l| , [37 



but without discussing the physical context of the problem. Our original motivation was to 
understand the apparent success of Tsallis statistics (i.e., the situations in which q > 1 or, 
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possibly also q < 1) in the realm of high energy collisions. It should be stressed that in this 
way we have addressed the interpretation of only very limited cases of applications of Tsallis 
statistics. They belong to the category in which the power laws physically appear as a conse- 
quence of some continuous spectra within appropriate integrals. It does not touch, however, 
a really hard case of applicability of Tsallis statistics, namely when zero Lyapunov exponents 



are involved ||38|| . Nevertheless, this allows us to interpret some nuclear collisions data in terms 
of fluctuations of the inverse temperature, providing thus an important hint to the origin of 
some systematics in the data, understanding of which is crucial in the search for a new state 



of matter, the Quark Gluon Plasma 0, 30]. 



5 Final remarks 

The are also other imprints of nonextensivity which we shall only mention. One is connected 



with recent analysis |12] of the equilibrium distribution of heavy quarks in Fokker-Planck dy- 
namics. It was demonstrated that thermalization of charmed quarks in a QGP proceeding via 
collisions with light quarks and gluons results in a spectral shape which can be described only 



by the Tsallis distrubution [09[] . On the other hand in |L3|] the quantum scattering processes 



(such as ttN and ttA) scatterings were analysed using Tsallis-like entropies and strong evidence 
for the nonextensivity were found there when analysing the experimental data on the respective 
phase shifts. On the boundary of really high energy collisions is the very recent application 
of the nonextensive statistics to the nuclear multifragmentation processes [4T5]. The other ex- 



amples do not refer to Tsallis thermostatistics directly, nevertheless it can be demonstrated 
that they are, at least approximately, connected to it. We would like to refer here to a recent 
attempt to study, by using the formalism of quantum groups, the so called Bose-Einstein cor- 



relations between identical particles observed in multiparticle reactions [O] and also works on 



intermittency and multiparticle distributions using the so called Levy stable distributions [42 



They belong, in some sense, to the domain of nonextensivity because, as was shown in [43 



there is close correspondence between the deformation parameter of quantum groups used in 



41 1 and the nonextensivity parameter q of Tsallis statistics and there is also connection between 



Tsallis statistics and Levy stable distributions Some traces of the possible nonextensive 



evolution of cascade type hadronization processes were also searched for in p5| . The quan- 
tum group approach |4l], [43! cou ld probably be a useful tool when studying delicate problem 



of interplay between QGP and hadrons produced from it. It is plausible that description in 
terms of g-deformed bosons (or the use of some kind of interpolating statistics) would lead to 
more general results than the simple use of nonextensive mean occupation numbers < n > q 
discussed above (for which the only known practical description is limited to small deviations 
from nonextensivity only). 

To the extend to which self-organized criticality (SOC) is connected with nonextensivity [T]] 
one should also mention here a very innovative (from the point of view of high energy collision) 
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application of the concept of SOC to such processes 



46 



To summarize, it has been demonstrated that multiparticle processes bear also some signs 
of nonextensivity observed in other branches of physics, which shows up only as small devia- 
tions from the expected behaviour. These deviations were already explained by invoking some 
additional mechanisms and, because of this, the use of g-statistics is not so popular or known 
in this field as in others discussed in 0. The advantage of the use of g-statistics is probably 
best seen from the information theoretical point of view. The new parameter q can be regarded 
then as a kind of compactification of all processes responsible for the actual nonextensivity 
into one single number f4~7| . This is also the point of view expressed in |48] where the new 



approach to quantum field theory based on Lorentzian, instead of Gaussian, path integrals has 
been proposed. It would allow to account for the possible deviations of pure stochasticity in a 
similarly most economical way when one introduces a single new parameter. This is, however, 
so far unexplored domain of research. 
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Figure Captions 



Fig. 1 Depth distribution of the starting points, dN(T)/dT, of cascades in Pamir lead chamber. 



Notice the non-exponential behaviour of data points (for their origin cf. |14j]) which can 
be fitted by Tsallis distribution @ with q = 1.3. (This figure is reproduced from Fig. 1 

of 01)- 

Fig. 2 The results for p T distribution dN{p T )/dpT- notice that q = 1.015 results describes also 
the tail of distribution not fitted by the conventional exponent (i.e., q — 1). This figure 
is reproduced from Fig. 3 of . 



Fig. 3 The example (Fig. 1 of ||) of transverse momentum distributions for e + e~ annihilation 
processes for differnt energies (shown in Figure). The inset shows region of small values 
of pt- The dotted line shows curve of constant T and q = 1. Other curves were obtained 
by fitting corresponding data with essentially fixed T ~ 110 MeV and q growing fast 
with energy to stabilize at ~ 90 GeV at value q = 1.2 (cf. || for details; this figure is 
reproduced from Fig. 1 of [|J). 

Fig. 4 The examples of the most probable rapidity distributions obtained by extending analysis 
of |T8|] (eq. (0)) to the nonexponential (q 1) distributions given by eq. flT3|) . The object 
(fireball, string,...) of mass M = 100 GeV decays into N secondaries of (transverse) mass 
m T = 0.4 GeV each. Figs, (a) and (b) show results for N leading to Feynman scaling 
{(3 = 0) or Feynman q-scaling (/3 g =i.3 = 0). Fig. (c) shows example of such iV that all 
P > 0. 



Fig. 5 Example of the nonextensivity in fluctuations: Fig. 1 of [[10] showing $ - measure of 
the kaon multiplicity fluctuations (in the tt~K~ system of particles) as a function of 
temperature for three values of the pion chemical potential. The kaon chemical potential 
vanishes. The resonances are neglected, (a) - results of |25| (in linear scale); (b) - our 
results for q = 1.015. 
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